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On ! Abstract 

Complete two-loop QCD corrections to b — ► c transitions are presented in 
the limit of zero recoil. Vector and axial-vector coefficients T)a V are calcu- 
lated analytically in the limit of equal beauty and charm masses, and a series 
appoximation is obtained for the general mass case. rjA is crucial for the 
determination of the absolute value of the Cabibbo-Kobayashi-Maskawa ma- 
trix element V c b- The two- loop effects enhance the one- loop corrections by 
\Q ■ 21%, removing a major theoretical uncertainty in the value of \V c b\- Including 

two-loop QCD effects and previously neglected electroweak and Coulomb cor- 
rections we find \V cb \ = 0.0379 ± 0.0020(stat) ± 0.0025(syst) ± 0.0011 (theory). 

\o 

ON 



Elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix are fundamental input pa- 
rameters of the Standard Model. Their precise measurements have been the subject of vast 
experimental efforts and will remain prominent issues in forthcoming projects, most notably 
the B-factories. The values of CKM matrix elements determine the sides of the unitarity 
triangle and their precise knowledge is essential for the understanding of the origin of CP 
violation, a major puzzle of the Standard Model. 

One of the directly measurable CKM parameters is the absolute value of Vcb- The exper- 
imental value can be extracted from decays of B mesons produced on the T(4) resonance 
(ARGUS Collab. 0, CLEO Collab. [§) or on the Z resonance (ALEPH CoUab. [§, 
DELPHI Collab. @J). 

\V c b\ can be obtained either from the total width of semileptonic B decays or from 
the zero-recoil extrapolation of the exclusive decay spectrum of B — > D*W, where I is an 
electron or muon (see || for a recent review). The merits of both methods and theoretical 
uncertainties have been discussed in ref. ||. The inclusive approach has the advantage of 
larger experimental statistics; the inherent theoretical error is mainly due to inaccurate 
knowledge of the quark masses which enter the decay width formula. This theoretical 
uncertainty already dominates the experimental error and it is not obvious that it can 
be significantly improved (see, however, a discussion in ref. @ and references therein). 
The exclusive method, on the other hand, benefits from recent advances |5],[S| in the heavy 



quark effective theory (HQET) It has been used to obtain the latest experimental 

result 



1 



\Vcb\ = 0.0385 ± 0.0021(stat) ± 0.0025(syst) ± 0.0017(theory) (ref. [4]). (1) 

The exclusive method can be summarized as follows: The recoil spectrum of the B meson 
decay is 

aw 

x(l + «.„^), (2) 

where u> is the product of the four- velocities of the B and D* mesons, and / is a known (see 
e.g. [§]) function of masses of observable particles (rather than of quark masses). HQET 
offers a model-independent value of the hadronic matrix element for the decay B — > D*lv at 
zero recoil, TiV)-, up to perturbative corrections, to be subsequently discussed. This point 
is not directly accessible in the experiment due to the vanishing phase space. Fortunately, 
|Vcft| 2 .F 2 (l) can be deduced by extrapolating the measured values at non-zero recoil, and, 
given the theoretical prediction for TiX), the value of |V^| can be obtained. The last factor 
in eq. (0) approximates the electroweak corrections In case of neutral B decays an 

additional factor (1 + 7ra) should be included in eq. @ (see |L4j] and references therein.) 
It represents a 2.3% enhancement of the rate due to the final state interaction between the 
lepton and the charged D* meson. 

In contrast to the decay into a pseudoscalar D meson, the prediction of HQET for 
B — > D* transition is free from 1/rrib^ corrections [15| by virtue of Luke's theorem |i~6| . The 
form factor J 7 (I) can be written as 

.F(l)=7 M (l + 5 1/m2 ), (3) 

where t]a = 1 + 0(a s ) describes the perturbative QCD corrections. The mass corrections 
Si/ m 2 of order 1/rriQ have been examined [PJTTfl. They are estimated [|l7j] to decrease the form 



factor by (5.5 ±2.5)% and their error is responsible for approximately half of the theoretical 
uncertainty in the value of \V c b\ quoted in eq. ([!]). 

The large remaining theoretical uncertainty is due to the unknown perturbative QCD 
corrections. The one-loop corrections are formally identical to those calculated in the context 
of muon decay |0| . For the heavy quark decays they give |l9| , |20lJl2 



+ (ff C F ,tf + 0(al), (4) 

where the dimensionless variable 5 = 1— rric/mb denotes the relative difference of quark 
masses, and rj\ represents the two-loop QCD corrections. The latter have been subject 
of vigorous controversy over the last few years. In the absence of an exact calculation, a 
renormalization group analysis has been performed [^TJ but its validity has been questioned 



in view of the small size of the logarithm of the mass ratio mb/m c [p2] , p3|| . The need for a 
full two-loop calculation of rf^ has been emphasized by many authors pl|-p4]| . The purpose 
of this paper is to provide this correction. 
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A calculation of QCD (or even QED) two-loop corrections to a fermion decay is in general 
very difficult. A full calculation has never been done, neither for the muon nor for a quark. 
However, it is at present possible to perform such an analysis at least at the zero recoil point. 
The advantage of this particular kinematical point is twofold. First, because of the phase- 
space suppression, there is no room for real radiation. Therefore, the virtual corrections are 
infrared finite. Second, the four vectors of the decaying and final quarks are parallel. 

These two features of the zero recoil configuration allow an exact analytic solution in the 
case of equal masses mj and m c ; in the general mass case one can construct an approximate 
solution in the form of a power series in the relative mass difference 5. In addition, the 
solution has a very useful symmetry with respect to the exchange mj m c . This non- 
trivial symmetry, valid only at zero recoil, helps to extract maximum information from the 
approximating series by accelerating its convergence. 

The two-loop QCD diagrams relevant to this calculation are shown in fig. 1. The Lorentz 
structure of the decay vertex is 

r M = 7^(1-75), (5) 

and the vector and axial-vector parts of this coupling are modified in different ways by the 
QCD corrections; at zero recoil they are parametrized by two functions, t]v,a 

In -> VvJn, 7m7s -> ?M7m75- (6) 
Although for the decay B — > D*W only the axial part is relevant, for the sake of completeness 



we also compute rjy. The one-loop result for rjy is known |T8|- pO| , p!2 

-5 



Vv 



71 




+ 



7T / 



(7) 



The one-loop QCD corrections contain only one color structure, proportional to Cp. At the 
two-loop level it is convenient to divide up the functions rj^ v into four parts proportional 
to various SU(3) factors (an overall factor Cp has been factored out): 



(2) 

Va'v 



CfVa,v + (Ca - ICp^fy 



(8) 



+T R N L ri L Ay + T R ril v . 

For a general SU(N) group C A — N; Cp — (N 2 - 1)/2JV; T R = 1/2. N L denotes the 
number of the light quark flavors whose masses can be neglected. The last term contains 
contributions of the massive quark loops, with b and c quarks. We neglect the top quark; 
its impact is suppressed by a factor ~ m 2 /m 2 . 

Among the eight coefficient functions in eq. (^) 7]\ v are already known EE \ . They 
corresponds to our diagram (f) in fig. 1 with a massless fermion in the loop. In the MS 
scheme (with /i = y/m b m c ), adopted also in the present work, they read 



V L A 



rib 



5 

24 
1 

24 



2-5 
5 

2-5 



ln(l 



5)+ 15 



ln(l -5) + 2 



(9) 
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The remaining six functions can be calculated exactly in the case of equal masses m& and 
m c . In this limit the momenta of the leptons in the final state vanish and the vertex function 
becomes a two-point function with a zero momentum insertion. Such propagator-like on- 
shell functions are known; a systematic method of their evaluation has been worked out in 
ref. p6|j27j| ; the underlying idea is the integration by parts method ||28| . This method has 



greatly simplified the two- loop QED calculation of g — 2 [^,0 . In ref. [29] the computation 
of two-loop functions with a low number of zero momentum insertions has been automated. 
For the purpose of the present calculation a new implementation of the recurrence algorithm 
3TJ was necessary; this is because of the necessity of computing two-loop functions with large 
number of zero momentum insertions. 

In order to go beyond the nib = m c limit we use the variable S as an expansion parameter. 
In the real world nib and m c are far from being equal; for the purpose of this work we take 
rrib = 4.8 GeV and m c = 1.44 GeV which yields 5 = 0.7. Coefficients of the expansion of 
r]Ay m <5 are two-point on-shell functions which can be computed using recurrence relations. 
In order to ensure good numerical accuracy we have computed ten terms in the S expansion 
for all diagrams in fig. 1. The analytic computation of the resulting integrals was feasible 



only thanks to the latest achievements in symbolic manipulation programs [32fl. 

The results we obtained are symmetric with respect to the exchange nib <-» m c , or 
6 — > —5/(1 — 5). The resulting fact that every other term in the 5 expansion can be 
obtained from the earlier terms provides a strong consistency check of our procedures. On 
the other hand, it is possible to rewrite the series expansion in a manifestly symmetric form. 
For this purpose we introduce a variable invariant with respect to nib *-> nic, P = 5 2 '/ (1 — 5). 
In terms of this variable our results can be written as a rapidly convergent series and we can 
safely neglect terms 0(p 5 ). Quadratic and linear logarithmic terms, implicitly present in 
this expansion P^|- p5|f2T[ ] , do not spoil the convergence because ln(m b /m c ) is of order unity. 



(2) 

For the axial- vector function ry A ' we find 



AF 143 1 o 1 2l 1 > , , 

rii F = vr 2 + -vr 2 n2 - -C 3 

IA 144 12 6 4 SV ; 

+p + vr 2 vr 2 ln2 + — C(3) 

H V576 1728 16 32 s v ; 

„ / 2509 121 2 \ , / 43 67 2 

+P 1 vr 2 + 7r 2 

H V 17280 8640 J H V 22680 967680 

4 / 17933 143 2 



+P" + 7T 

H V 50803200 11612160 
F 373 1 2 

/377 1 2 ^ : ^ 2 ( 29 1 2 

'648 + 360 7r 
o / 227 1 2 \ , ( 649 1 2 

+P 7T 2 + p 4 1 7T 2 

H V37800 2520 J H \ 672000 15120 
h 115 2 2 



+P 1,576 - 24^ ) + " 



/529 107 2 \ 2 ( 337 137 2 



72 144 71 ) +P ri44 + 576 W 
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, i 255313 197 2 \ 4 / 1957573 1 2 

-p 3 1 tt 2 + p 4 1 7T 2 

; 1 75600 576 / H V 3175200 16 



For the corrections to the vector current we find: 

^-"(I-^ 2 -^ 2ln2+ l« 3 ) 

2/107 1 2 \ , 3 / 41 31 2 



(10) 



+P 1 7T + p° IT 

1 1 5760 288 / H V 10080 46080 
, , 13927 169 2 



+ P" + 7T 

' 1 16934400 1290240 

P /553 5 2 \ 2/ 227 1 2 
= p ^ - — vr 2 + p 2 -— - + — vr 2 



„3 



V576 72 / r V 4320 360 
/ 251 1 2 \ 4 / 7537 1 



+p vr + P 1 7i" 

H V33600 2520 J H \ 6048000 15120 

H /197 13 2 \ 2 / 701 19 2 



^ ='072-48* J +P 0720 + 192" 

, / 2851 55 2 \ .( 93227 1 a 

+p 1 vr 2 + p 4 1 7T 2 

r V 1008 192 J p V 151200 16 



(11) 



We note that all QCD contributions to rjv vanish at m& = m c (p = 0), in consequence of 
vector current conservation. 

So far we have not discussed the renormalization procedure which led to the results in 
eq. (p"0|, 11). For the external quark legs we used the two- loop quark wave function renor- 



malization constant computed in [[3(J. Vanishing of the terms independent of p in eq. (JTTJ) 
serves as an independent check of the complicated expressions given in [j36| |. The diagrams 



(bl) and (b3) require mass counterterms. For these we adopted the on-shell condition. Our 
results are therefore in terms of the pole masses mj, and m c . For the coupling constant renor- 
malization we used the MS condition with the renormalization scale at the geometric mean 
mass fi = ^fm^m~ c . It must be noted that the symmetry m& «-> m c is in general valid only 
for the unrenormalized diagrams. If the coupling constant were normalized at a different \i 
the final result (0,0) would not be symmetric. The geometric mean point is the only one 
preserving that symmetry. 

Numerically, the two-loop corrections evaluate to 

rjf = -0.586(2) x C F - 0.909(2) x (C A - 2C F ) 
+0.145 x T R N L - 0.155(4) x T R = -0.944(5), 

v p = 0.395(2) x C F - 0.168(2) x (C A - 2C F ) 

-0.010 x T R N L + 0.107(2) xT R = 0.509(5). (12) 

It is interesting to compare these results with an estimate based on the subset of corrections 



of order a 2 /3 with /3 = 11 — ^nj. With n/ = 4 one gets |25 

rjW = _o.908, r$ = 0.061. (13) 
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In the axial-vector case the agreement with the full two-loop calculation is quite good. The 
estimate fails badly in the vector case; this is probably because of an accidental numerical 
cancellation in (|13|) which makes the estimate of rjy^ very small. The full one and two- loop 
corrections tend to give corrections to r] V with approximately half the magnitude of those 

to 7]A- 

Adopting a s ( v /mf,m c ) = 0.24, our full two loop calculation leads to the total values of 

Va,v 

Va = 1 - 0.033 - 0.007 = 0.960 ± 0.007, 

Vv = 1 + 0.018 + 0.004 = 1.022 ± 0.004. (14) 

Since the perturbative series in QCD is asymptotic, the uncertainty in the values of t]a,v 
has been estimated by the size of the last computed terms. The central value we obtain for 
rjA is consistent with the value given by Neubert, ^(ref. ||) = 0.965 ± 0.020, which was 
adopted in the recent experimental studies. Our result reduces the error bar by a factor 3 
and removes a major source of the theoretical uncertainty in \V c b\. This error can perhaps be 
further decreased by choosing a different renormalization scheme, e.g. the ^/-scheme. This 
possibility will be examined in a future work. However, with our estimate of the perturbative 
two- loop corrections to tja, the uncertainty in the zero recoil form factor is dominated by 
the error in the l/rriQ corrections; we adopt here the value 5y m 2 = —(5.5 ± 2.5)% fll7|j8|| 
(for a recent discussion of these corrections see f3~?|| .) Putting this result together with our 
rjA = 0.960 ± 0.007 we find for the zero recoil form factor 

= ri A (l + (J 1/m a) = 0.907 ± 0.026 (15) 

We use the latest experimental data from the recent DELPHI analysis for the B° — > 
D*~l + v decay rate. We include the electroweak correction and the Coulomb enhancement, 
as discussed after eq. ||. They enhance the rate by 1.3% and 2.3%, respectively. Altogether, 
we find 

\V cb \ = 0.0379 ± 0.0020(stat) 

±0.0025(syst) ± 0.0011 (theory). (16) 

We see that with the decreased theoretical uncertainty further improvement in statistical 
and systematic accuracy can significantly increase the precision of \V c b\ and bring us closer 
to overconstraining the unitarity triangle. 
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FIGURES 




(e) (f) 

FIG. 1. Two-loop QCD corrections to the b — ► c transitions at zero recoil. Symbols ® mark 
the coupling of the virtual W boson to the quark line. 
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